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FIG. 1: The spectral region containing the 3A S molecular 
vibration, Raman spectra (solid lines) and Infrared spectra 
with polarization parallel to the crystallographic 6-axis (dot- 
ted lines), for a selection of pressures at room temperature. 
The spectra are offset for clarity and the vertical lines are a 
guide to the eye. Note the substantial softening effect of the 
increased electron-phonon interaction on the IR spectra with 
respect to the Raman spectrum. 

There are two approaches to modelling the electron- 
molecular vibration coupling interaction and associ- 
ated mode softening in organic charge transfer sys- 
tems, the "phase-phonon" theory and the "dimer charge- 
oscillation" model. The phase-phonon theory approaches 
the problem from an electron band approximation where 
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the charge carriers are naturally delocalized. The al- 
ternative dimer charge-oscillation theory is formulated 
in real space such that the electronic behaviour is im- 
plicitly localized. In reality the electronic properties of 
k-(BEDT-TTF) 2 Cu(SCN) 2 are somewhere hetween the 
extremes of a localized and itinerant systemic. For this 
reason we apply both the dimer charge-oscillation and 
phase-phonon models to analyze the data. 

This paper is organized as follows. In Section |l| 
we describe the vibrational spectra of k-(BEDT- 
TTF) 2 Cu(SCN) 2 as observed by both non-resonant Ra- 
man scattering and polarized IR reflectivity. Section III 



describes the methods of analysis, specifically how the 
pressure dependence of the electron-phonon coupling is 
extr acted. The first subsection of the Discussion, Sec- 
tion IV A, contains a comparison of the quasiparticle 



masses that result from the differ ent methods of analy- 
sis. Subsequently, in Section IV B, in light of the pressure 
dependence of the electron-phonon coupling, the impli- 
cations for the superconducting pairing mechanism are 
discussed. Conclusions are summarized in Section |y|. 
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FIG. 2: The room temperature diamond/sample IR re- 
flectance offset by pressure for polarization parallel to the b- 
axis. Measurement (solid lines) and phase-phonon fit (dashed 
lines). The inset magnifies the spectral region containing the 
C=C mode. 



II. DESCRIPTION OF THE OBSERVED 
VIBRATIONAL MODES. 

Four modes are observed in both the high-pressure 
room-temperature Raman and IR spectra (see Figures ^ 
and |^ and Table |). They are labeled with subscripts in- 
dicating the atoms/bond predominantly involved in the 
vibrationEa. In order of increasing frequency they are 
the C-S mode originating from the BEDT-TTF 60B 3g 
asymmetric vibration, the central C=C mode originat- 
ing from the BEDT-TTF 3A g symmetric vibration and 
two Cu(SCN) 2 anion modes. Determining the pressure 
dependence of the central C=C mode in the IR spec- 
trum is not straightforward due to the overlap of several 
modes with— varying pressure dependences. In a recent 
publications, IR spectra from the deuterated salt have 
been used to model this Fermi resonance and extract the 
linear pressure dependence of the overlapping modes (the 
central C=C mode and two C-H modes). The IR and Ra- 
man frequencies and the first order pressure dependence 
of all modes discussed here are given in Table Q. 

The Cu(SCN) 2 anion modes do not exhibit any evi- 
dence of electron-phonon coupling, i. e. the difference be- 
tween the infrared and Raman frequencies for the anion 
modes is accounted for by purely vibrational couplingEI. 
The zone center frequency separation, Ao> s , between the 
symmetric (Raman active) and antisymmetric (infrared 
active) combinations of the molecular vibrations is de- 
termined by the. frequency of the optical branch of the 
lattice modesEO. The larger observed Aw s for the anion 
mode in the 6-axis response, for which the lattice mode 
is stiffer (see Table |), is a further confirmation of the 
lattice mode assignment given inQ. 
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FIG. 3: The room temperature diamond/sample IR re- 
flectance offset by pressure for polarization parallel to the c- 
axis. Measurement (solid lines) and phase-phonon fit (dashed 
fines). The inset magnifies the spectral region containing the 
C=C mode. 

III. ANALYSIS. 

Electron-molecular vibration coupling can be charac- 
terized by a set of linear coupling constants, g a , which pa- 
rameterize the interaction between the a th normal mode, 
of frequency uj a , and the molecular orbital in which the 
radical electron or hole residesEj. For a donor molecule 
g a is the rate of change of the energy of the highest occu- 
pied molecular pcbital with respect to the normal mode 
coordinate, Q Q ,li3 
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Mode Assignment Raman Infrared 6-axis Infrared c-axis 
]l| cm' 1 +(%GPa~ 1 ) cm' 1 +(%GPa~ 1 ) cm' 1 +(%GPa~ 1 ) 



wot 








2910 


+ 4.0 


2390 


+ 4.0 


wcs 


B 3s 


886.2 


+ 0.85 


883.5 


+ 0.71 


873.6 


+ 1.0 


wcc 




1467.7 


+ 0.4 


1290 


+ 2.5 


1210 


+ 5.5 


WCH! 








1181 


+ 0.5 


1177 


+ 0.5 


WCH 2 


Bi u 






1290 


+ 0.5 


1281 


+ 0.5 


^anionl 


CN 


2064.6 


+ 0.1 


2067.4 


+ 0.1 


2065.6 


+ 0.15 


^anion2 


CN 


2106.3 


+ 0.2 






2109.3 


+ 0.2 



TABLE I: Raman and IR frequencies and pressure shifts taken from! 



To obtain these coupling constants experimentally, an 
appropriate model has to be applied to the data. 
Two alternative approaches exist in the literature: 
the "phase-phonon" theory and the "dimer charge- 
oscillati on" .model. Both theories were initially developed 
by RiceMliZl to model the infrared spectrum of semicon- 
ducting organic charge-transfer systems. However, they 
have also been applied to a range of reduced dimen- 
sionalit y organic , charge , transfer salts, including metallic 
systemsMliifcaM'Efi 

The phase-phonon theory approaches the problem 
from an electron band approximation where the charge 
carriers are naturally delocalized. A potential disadvan- 
tage of this approach is that it is a one-electron theory, 
i.e. all electron-electron correlation effects have. to be 
included within an effective mass approximationEJ. 

The alternative dimer charge-oscillation theory is for- 
mulated in real space but is limited to a finite number of 
molecules. The infrared response of the system is calcu- 
lated from a superposition of the isolated units with no 
significant charge transfer between the units. For /c-phase 
BEDT-TTF salts the approximation is that the charge 
carriers are predominantly localized on the dimers which 
are weakly interacting with each other. Although this im- 
plies an insulating system, with the response of delocal- 
ized charge carriers parameterized separately, electron- 
electron correlations such as the onsite Coulomb repul- 
sion can be included explicitly^. 



A. The dimer charge-oscillation model. 

This subsection outlines_lhe derivation of the dimer 
charge-oscillation modelt3'E£l. Not only does this aid dis- 
cussion relating to the validity of applying this model to 
k-(BEDT-TTF) 2 Cu(SCN) 2 , but it also lends theoretical 
weight to the concept of electron coupling to the anti- 
symmetric combinations of molecular vibrations and not 
to the symmetric combinations. This concept is the key 
to comparing Raman and infrared data in order to probe 
the electron-phonon coupling. 

Each dimer unit is modelled by the following Hamilto- 



nian (H = 1)00, 

H = H c + H v + Y^ njgaQcj. (2) 

The first two terms describe, respectively, the electronic 
states and the molecular vibrations in the absence of vi- 
brational coupling; j labels the molecule and a the nor- 
mal mode of vibration. Rice's formulation includes only 
the totally symmetric (A g ) molecular modes. However, 
later theories have extended the theory to internal modes 
of other symmetriescj. 

The electronic term, H e , is of the form, 

3 = 1:2 

where nj = c ] <r c 3,v ^ s ^ ne occupation number opera- 
tor for the molecular orbital of energy Eq. c, and Cj i<7 
are the Fermion creation and annihilation operators for 
an electron of spin oEj'EI The exact form of the potential 
term, V , is arbitrary. 

The phonon term, H u , is of the co nven tional form for 
a harmonic oscillator of frequency wJISEJ'Ej, 

ff, = £w a (&5,Aa + ~)> (4) 

j-a 

where b\ „ and b,- a are the Boson creation and annihila- 
tion operators. Their sum is the displacement operator 
for the a th normal, mode and their difference is its mo- 
mentum operatorc3, i.e. 

Qj,a OC (bj a + bj ta ). (5) 

Because the dimer has been chosen as the base unit for 
this system it is appropriate to introduce dimer-mode 
coordinates 

s a = -j={Qi, a + Qi,a) (6) 

and 

(lot = ~J=(Ql,a - Q%,ct) (7) 
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which are symmetric and antisymmetric combinations of 
the molecular-mode coordinates respectively. 

The third term in the Hamiltonian is a linear pertur- 
bation arising from the electron-molecular vibration cou- 
pling. Commutation of the dimer-mode coordinates with 
the system Hamiltonian yi elds .two equations of motion 
for the a th molecular moddI3lri] 



(s a ) + Lol(s a ) = -V2g a u a (ni + n 2 ) 



and 



(8) 



(9) 



The right-hand side of the equation of motion for the 
symmetric combination of molecular modes is depen- 
dent on the total charge density, (iVtot) = ("l + "2), 
which for an isolated dimer is constant, i.e. the sym- 
metric combination of modes is uncoupled from the elec- 
tron system. Their unperturbed frequencies, u> a , are 
therefore obtained experimentally from Raman scatter- 
ing experimentsEj. 

The right-hand side of the equation of motion for the 
antisymmetric combination of molecular modes is de- 
pendent upon the electric dipole moment of the dimer, 
(p) = ^ea(ni — n 2 ), where e is the electron charge and 
a is a vector perpendicular to the interface between the 
two molecules, whose magnitude is of the order of the 
intermolecular separation. Therefore the antisymmetric 
combination of modes is directly coupled tcuthe electron 
dipole moment and is hence infrared activeEj. For a sim- 
ple two-state electron Hamiltonian these modes couple 
to the charge-transfer excitation between the ground and 
excited states, |1) — > |2). 

Because the interaction part of the Hamiltonian is of 
the form 



(10) 



where ip is a scalar potential, the expectation value of 
the operator p which characterizes the response of the 
system to the perturbation may be calculated using lin- 
ear response theoryt^Oo. The ground state expectation 
value of the Fourier components of the dipole moment, 
(p(u>)), and hence the complex conductivity, <j(u>), of the 
dimer charge-osci llati on system are given by the linear- 
response functionl!3E3, 



cr(u>) = —iu> 



e 2 a 2 / 1 



E 



2V \ XM „ w q - ^ 2 - ' luJ la 



(11) 

Here V is the molecular volume and x( w ) t ne reduced 
charge transfer polarizability; 



xH = I(2|p|i 



1 2 2ucT 
UIq t — OJ^ — iwl 



(12) 



where (2|p|l) is the matrix element for the charge- 
transfer transition of frequency wct- The electron and 



phonon damping factors, T and j a , are introduced phe- 
nomenologically to account for the observed infrared 
linewidths. It should be noted that this perturbation 
will alter the frequency at which the coupled molecu- 
lar vibrations are observed in the infrared spectrum; for 
uj a > wct their frequencies will be red shiftecQ 

The crux of this subsection is to address how the 
dimer charge-oscillation model may be used to extract 
the electron-phonon coupling constants by means of ex- 
periment. Several approaches have been proposed to re- 
duce the possible parameter space encountered when fit- 
ting infrared reflectance data with this model. It has been 
notediS that the inverse of the real part of the conductiv- 
ity, calculated using equation [ll], consists of a constant 
background arising from the charge-transfer band with 
peaks occurring at the unperturbed frequencies of the 
vibrational modes, u) a . In principal, if this model accu- 
rately describes the data, Kramers-Kronig analysis could 

be used to calculate Re f^^yj and thus determine the 

bare-mode frequencies. However, because of uncertain- 
ties in the Kramers-Kronig procedures Raman measure- 
ment of the unperturbed frequencies are more reliable. 

Analysis of the poles and zeros of the response function 
(Equation [ll]) leads to the following approximate relation 
between the observed frequency of the infrared phonon 
bands, Q a , their unperturbed frequencies, oj a , the fre- 
quency of the coupling charge-transfer band, loct, and 
the dimcnsionless electron molecular vibration coupling 
constants, A Q li3'E3; 



— x a - 



CT 



(13) 



This expression is valid if the bare frequencies of the 
modes are well separated from each other or X a suffi- 
ciently small so that each mode may be treated individ- 
ually. For k-(BEDT-TTF) 2 Cu(SCN) 2 the two observed 
modes with a finite A Q (the C=C and C-S modes) are 
well separated from each other. Note that the C=C mode 
overlaps with two C-H modes as discussed in Section H 
however theJack of spectral weight associated with the 
C-H modesEj indicates that their electron-phonon cou- 
pling is negligable in this case. Within the dimer charge 
oscillation model X a is related to g a viaE3 



2gl 



(14) 



B. The electron-phonon coupling constant and its 
pressure derivative derived from the dimer 
charge-oscillation model. 

Equation |l^ is used to calculate an electron-phonon 
coupling constant, A Q , for each mode from its Raman 
and infrared frequencies, u a and Q , given knowledge a£ 
the energy, ojqt, of the coupling charge-transfer bandl3't£l 
(see Table ||) . The dominant source of error in this cal- 
culation is due to the width of the charge-transfer band. 
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There also exists the possibility of coupling to two differ- 
ent charge-transfer transitions. The C=C mode, being 
an antisymmetric combination of A g molecular modes, 
couples to .the intradimer charge transferee. It has been 
calculated,E3 however, that the antiphase combination of 
T$3 g molecular modes couples to charge transfer perpen- 
dicular to the intradimer direction. This suggests that 
the C-S mode couples to transitions between the lower 
and upper branches of the same band, not interband 
electronic transitions. However, it should be noted that 
the electron-phonon coupling constants for the C-S mode 
have been calculated using the same ujqt as the C=C 
mode because it is impossible to distinguish the contribu- 
tions to the infrared spectrum from different transitions. 



The pressure derivative of Equation 13 



din A a 
dP 



= 2- 



dine 



dP 



d in Q a 
dP 



,2 



d In t^cT d in to a 



dP 



dP 



(15) 



has a weaker dependence on the value of loqt- How- 
ever the dominant source of error in Equation [l^ arises 
from the difficulty in determining an accurate pressure 
dependence of the broad charge-transfer band. Fitting 
the infrared reflectance data with a simple Drude-Lorentz 
model! provides an estimate for the pressure dependance 
of us&r. Using a value of 4±4 %GPa _1 for both b- and o 
axesEl, results in an error in the pressure derivative of the 
coupling constant for the C=C mode of ±3.5 %GPa _1 
for the 6-axis and ±5 %GPa _1 for the c-axis. The er- 
ror in the pressure derivative of the coupling constant 
for the C-S mode is larger (8 %GPa _1 for the &-axis and 
28 %GPa _1 for the c-axis) because Acs is so small. 



Mode 


A„ 


d In Aj, 
dP 

^oGPa" 1 


A c 


d In A c 
dP 




0.01(1) 


47.9 


0.02(1) 


-11.1 




0.17(1) 


-11.9 


0.20(1) 


-17.3 


k-'anionl 














^anion2 













TABLE II: dimensionless electron-phonon coupling constants 
and their pressure derivatives calculated using the dimer 
charge-oscillation model. Note in this table the uj subscript 
referes to the mode and the A subscript to the polarization 
direction. 



Because the dimer charge-oscillation model is formu- 
lated for a localized system, A Q calculated in this man- 
ner parameterizes the strength of interaction between 
the phonons and the intradimer charge transfer. On the 
other hand the electron-phonon coupling constant asso- 
ciated with phonon-mediated superconductivity param- 
eterizes the strength of interaction between the phonons 
and the delocalized charge carriers. Therefore to use 



this model to draw any conclusions regarding the super- 
conducting mechanism requires the assumption that the 
electron-phonon coupling for the intradimer charge trans- 
fer will, to first order, have the same pressure dependence 
as the electron-phonon coupling for the interdimer charge 
transfer. This assumption is reasonable because it is the 
same molecular orbitals that are responsible for intra- 
and interdimer wave function overlap. At the T-point 
(the lone_wavelength limit) neighboring dimers vibrate 
in phaseEj and as a result, an antiphase combination 
of molecular vibrations that modulates the intradimer 
wavefunction overlap will also modulate the interdimer 
wavefunction overlap. 



C. The phase-phonon model. 

Phase-phonon theory essentially describes the same 
phenomenon as the dimer charge-oscillation theory (see 
equation |^). However, because it is formulated in re- 
ciprocal space the model is easily extended to metallic 
systems. This model is not just limited to the dimer 
unit, but attributes the infrared activity of the coupled 
phonon modes to phase oscillations of the spatial charge 
density, induced by the electron-phonon interaction!!!. 

The system Hamiltonian is similar to that used in 
the dimer charge-oscillation theory, and includes an elec- 
tronic term, H C1 a phonon term, H v , and a linear cou- 
pling term, H{. 



H = H C + H U + -^=Y^ 9 a Qa{q)p 



(16) 



a: , q 



The electronic terms has essentially the same form as (|^) , 
H c = J2Ekcl„c k ^ + V(p qo +p- qn ), (17) 

except that the operators create or annihilate a parti- 
cle in a fc-state, not at a real space location, i.e. it de- 
scribes a system of conduction electrons moving in a weak 
periodic potential, V, of wavevector q . The operator 
P q = J2k c l+ q c k creates an electronic-density fluctuation 
of wavevector q. 

The phonon term is identical to that used in the dimer- 
charge model except the phonon operators also act on 
fe-states. The electron-molecular vibration coupling is 
included as a linear perturbation, in this case linking 
the Fourier transform of the mode displacement vector, 
Qa{q) oc (b a (q) + &£,(—<?)), to an electronic density fluc- 
tuation of the same wavevector, q. 



Hi = 



77= 9 a c{ +q c k (b a (q) + bl(-q)). 

a,q,k 



(18) 



N is the density of conduction electrons. ■— ■ 

For the case initially considered by RiceJlH go = 2fci? 
(twice the Fermi wavevector) such that the periodic 
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potential induces the conduction electrons to condense 
into a charge-density-wave state. However, the deriva- 
tion that follows is equally applicable to a metallic sys- 
tem. The frequency-dependent conductivity for an in- 
sulating svstem is due simply to interband electronic 
transitionsEj, 



a(cj) = 7 Z[f(x)-f(0)\ 



(19) 



where the plasma frequency, u> p , parameterizes the band 
curvatures. /(0) = 1 and 



in 



s 

1+S 



2Sx 2 



where S 



and x 



2A 



(20) 



with Aq equal to 



the band gap. Including the electron-phonon coupling 
modulates the onsite electron energy and hence the band 
gap, so that A = Ao + ^ Q A Q e"^ Q , where A Q and <j) a 
are the amplitude and phase respectively of the distor- 
tion potentials which are proportional to g a . The sin- 
gle electron contributions to cr(uS) are of the same form 
as the uncoupled case but with x = In addition 

to the single-electron contributions there are collective 
contributions associated with oscillations in the phases, 
(j) a , of the combined lattice-charge distortions. Collective 
modes associated with oscillations in amplitude, A Q , also 
occur: however, they preserve dipole momeetpsuid hence 
do not contribute to the dielectric responseC-JEll. 

Including the collective mode contribution, the 
frequency-dependent conductivity is given bye-3 



a(u) = ~~~ [f(x] - /(0) - x 2 (f(x)) 2 \D^(x)] . (21) 

Drj, is a phonon-like propagator for the phase oscillations 
given by 



1 _ 1 



1 



^ + \x2f(x) 



(22) 



Dq is the same single-phonon propagator used in the 
dimer charge-oscillation model, 



D Q = - 



(23) 



In the phase-phonon model the dimensionless coupling 
constant, A is related to g a via 



Aa — 



Hkp)g 2 



(24) 



where N(fep) is the density of electron states at the Fermi- 
surface and-the total electron-phonon coupling constant, 

A = E a aH. 

It is necessary to include an electronic damping term, 
r, to account for the finite electron lifetime if the model 



is to be applied to a metallic system, i.e. when a large 
plasma frequency is required to account for the intraband 
processes. In this case f(x) is replaced by f(x + iT)l23. 

For uj a < 2A decay of the collective mode via excita- 
tion of a real electron-hole pair is energetically impossi- 
ble. As a result each mode contributes a sharp absorption 
band to the infrared spectrum whose width is determined 
solely by the lineshape of the uncoupled phonon mode. 
For oj a > 2A the collective modes become damped via 
electron-hole excitation, appearing as indentations in the 
electronic background^. 



D. The electron-phonon coupling constant and its 
pressure derivative derived from the phase-phonon 
model. 

The electron-phonon coupling constant that occurs in 
the phase-phonon model parameterizes the degree to 
which each molecular vibration modulates the energy of 
the electronic band structure, i.e. it parameterizes both 
the inter- and intraband electron-phonon coupling. 

To calculate A a for each mode the full diamond/sample 
IR reflectance spectrumQ is fitted at each pressure-point. 
The model function for this fit consists of the phase- 
phonon model plus a high frequency contribution to the 
dielectric constant, £oo, and a highly damped Lorentzian 
oscillator to account for the anomalous electronic damp- 
ing at room temper atureE3. The plasma frequency, elec- 
tronic damping and perturbed and unperturbed inter- 
band transition frequencies are treated as free parame- 
ters. As with a simple Drude response, the plasma fre- 
quency parameterizes the density of states and band cur- 
vature, but in this case the electronic damping controls 
the relative spectral weights of the inter- and intraband 
processes. 

Four phase-phonon collective modes (en = 1 — ► 4) are 
used to model the spectrum, one for the C-S mode and 
three for the C=C mode and its Fermi resonance with 
the C-H modes. It was found that the line shape of the 
Fermi resonance could be reproduced in an-.analogous 
manner to the Green's function model used inlla, i.e. with 
one strongly coupled C=C mode and two C-H modes 
with negligible electron-phonon coupling (see the inset of 
Figures g and |3|) . The pressure dependence of the modes' 
Raman frequencies and the pressure dependence of the 
dips from the Green's function model were used to fix the 
frequency of the collective modes. Their damping and 
coupling strengths were left as free parameters, however. 

In this model the degree of mode softening associated 
with the electron-phonon coupling is predominantly pa- 
rameterized by the perturbation to the interband transi- 
tion, i.e. the difference between the perturbed and un- 
perturbed transition frequencies. The coupling strength, 
Aq,, controls the mode softening to a lesser extent, mainly 
parameterizing the spectral weight associated with each 
mode. Because of this, correctly determining A relies on 
accurately modelling the background reflectance in the 



7 




FIG. 4: The pressure dependence of the free parameters of the phase-phonon model. Note that the linear fits indicate the first 
order pressure dependence used for the batch fitting procedure. Solid lines and symbols are used for polarization parallel to 
the fe-axis and dashed lines and hollow symbols for polarization parallel to the c-axis. (a) the plasma frequency and electronic 
damping, (b) the difference between the perturbed and unperturbed interband transition frequency, (c) the high frequency 
dielectric constant, (d) the strength of the Lorentzian oscillator, (e) the frequency and damping of the Lorentzian oscillator, 
(f) the damping of the C-S and C=C modes. 



spectral region containing the modes. To this end, it 
was found that including the heavily damped Lorentzian 
to help model this material's anomalous Drude response 
drastically reduced the scatter in the A data. The three 
parameters describing the Lorentzian, uj e ,j e and Ae e , 
were also free. 

The lineshape of the anion modes is most accurately re- 



produced using uncoupled Lorentzian oscillators. Owing 
to their lack of electron coupling and their limited con- 
tribution to the spectrum, parameters describing them 
were not included in this model. Figures ^ and g indi- 
cate the quality of the fits for a representative selection 
of pressures. 

Manually fitting 31 spectra, each consisting of 3300 
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FIG. 5: The pressure dependence of the electron-phonon coupling constant. Solid squares and lines for polarization parallel to 
the 6-axis and hollow circles and dashed lines for polarization parallel to the c-axis. (a) for the C-S mode and (b) for the C=C 
mode. 



Parameter 


b-axis 


value 


c-axis 


value 




cm" 1 


+%GPa" 1 


cm -1 


+%GPa" 1 




(6816±120) 


+ (13±2) 


(5190±66) 


+ (19±2) 


r 


(815±16) 


+(45±3) 


(781±12) 


+(45±3) 


A 


(1132±13) 


+(10±2) 


(743±7) 


+(17±2) 


A 


(1407±14) 


+(6±2) 


(1064±8) 


+(6±1) 


A- Ac, 


(275±3) 


-(7±2) 


(320±3) 


-(19±2) 


7cs 


(15±0.5) 


-(26±7) 


(13±0.5) 


+(8±8) 


7cc 


(68±2) 


-(21±5) 


(101±5) 


-(18±6) 


Ul e 


(1905±25) 


-(28±2) 


(1257±21) 


-(12±2) 


7e 


(1214±51) 


+ (2±6) 


(1065±38) 


-(21±2) 






+%GPa" 1 




+%GPa" 1 


Ae e 


(0.90±0.15) 


+(388±22) 


(1.00±0.10) 


+(272±15) 


£00 


(3.9±0.1) 


+(11±3) 


(3.4±0.1) 


+(18±3) 


Acs 


(3xl0 _3 ±lxl0 _4 ) 


-(32±5) 


(3xl(T 3 ±lxl(n 4 ) 


-(25±4) 


Acc 


(0.080±0.002) 


-(55±4) 


(0.073±0.002) 


-(44±3) 



TABLE III: Phase phonon model parameters and their first order pressure shifts. 



data points, with an expression containing 12 free pa- 
rameters, for each polarization, is both immensely time 
consuming and subject to systematic human error. To 
avoid this a batch fitting procedure was employed. It con- 
sists of accurately fitting a limited number of data sets 
to determine approximate trends for the free parameters. 
These are then used as the the starting conditions for an 
automatic least-squares fitting procedure of all data sets. 
Each fit is then checked visually for acceptance of its pa- 
rameters. 

The value of the free parameters and their first order 
pressure shifts extracted from Figures f| and [3] are con- 



tained in Table [II 



The difference in the values of the 
coupling constants, A, obtained from the phase-phonon 
model (TABLE III) and the dimer charge-oscillation 



model(TABLE ||) are predominantly due to the different 
values for the interband transition used in each calcula- 
tion. 



IV. DISCUSSION. 

As well as vibrational features the room tempera- 
ture IR spectrum of k-(BEDT-TTF) 2 Cu(SCN) 2 contains 
broad electronic features; a heavly damped Drude-like 
response in the far-infrared and a mid-infrared (MIR) 
'hump'. Regarding the MIR 'hump' the literature only 
agrees to the extent that it originates-fcom, excitation 
across a gap in the electronic spectrumE3't3'E3. The sim- 
plest picture is to assign it to a gap in the single electron 
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band structural. This model is improved by including 
electron-electron interaction effectsL-X In this case the 
'hump' originates from a combination of the single elec- 
tron band gap and the onsite correlation energy. Which 
is considered dominant, and which the perturbation, de- 
pends upon how correlated the system is believed to be. 
Analogously, inclusion of the electron-phonon interac- 
tion modifies the electronic spectrum. The phase-phonon 
model includes this explicitly as a relatively small per- 
turbation to the band gap. In the limit of large electron- 
phonon coupling |the strain-held is considered to localize 
the charge carrieiH, in which case the small polaron bind- 
ing energy will be the band gap dressed by the electron- 
phonon interaction. All these theories share the common 
assumption that IR active interband electronic processes 
are subject to remormalization by many-body effects. In 
the following subsection we do not attempt to identify 
the origin of the renormalizations, but present evidence 
of their influence on the pressure dependence of the elec- 
tronic contributions to the IR spectrum. 



A. The pressure dependence of the carrier effective 

mass. 

This subsection focuses on the pressure dependence of 
the carrier effective mass extracted from modeling with 
the phase-phonon theory. The plasma frequency, cj p , pa- 
rameterizes the charge carrier density, N, and effective 
mass, here denoted ml 



pp' 



Ne z 



eom 



(25) 



pp 



For our actual pressure range the band filling is indepen- 
dent of pressure and the carrier density only has a small 
pressure dependence arising from the reduction in unit 
cell parametersoO. 

For a plasma frequency dominated by intraband pro- 
cesses the pressure dependence of m* p arises solelyfrom 
the pressure dependence of the band parameterstm, i.e. 
the increase in band width and curvature as the wave 
function overlap increases. We argue that if the plasma 
frequency also parameterizes the interband processes it 
will contain contributions from quasiparticle ma ss ren or- 
malization. The first point to note (see Section [ill C|) is 
that the plasma frequency in the phase phonon model not 
only describes screening ability of the free carriers (the 
intraband response) but also the density of states and 
band curvatures associated with intraband absorption. 
Thus following the preceeding logic the mass m* p will to 
some extent be renormalized by many-body effectsO. 

The second point to note is that using the pressure 
dependence of the Raman modes to constrain the mode 
frequencies for the phase-phonon model significantly re- 
duces the parameter space encountered when modeling 
the IR response of this system. This is not only because 
it eliminates the mode frequencies as free parameters, but 



because the softening of the modes is dependent upon the 
perturbation to the band gap, A — Ao. In this manner 
the mode frequencies provide further constraint to the 
parameters describing the electronic response. The pa- 
rameters describing the electronic response are thus be- 
lieved to have greater validity than those obtained from 
a simple Drude-Lorentz ffo3. 

The pressure dependence of the effective masses m* p 
obtained using the phase-phonon theory is shown in Fig- 
ures [| and 0. For both polarization directions m* p 




0.5 1.0 1.5 

Pressure (GPa) 



2.0 



FIG. 6: 



The mass m,p P 



for light polarized parallel to the b- 



axis, deduced from the phase-phonon fit to the optical data. 

For P < 0.5PGa, dm**/dP ' 2.2m. 

-80%/GPa. For P > 



b /dP ~ - 
0.5PGa, 



and dlnm* p 



d\nm* p b /dP ~ -5%/GPa. 



dm* p b /dP 



b /dP ~ 
0.9m e and 



decreases with the application of pressure. For both po- 
larizations this decrease is clearly nonlinear and this is 
highlighted by linear fits to the low (< 0.5 GPa) and 
high pressure regions (> 0.5 GPa)(see Fig. ^ and [?]). 
These pressure ranges were chosen because the 6-axis 
m pp b clearly exhibits a change of slope at 0.5 GPa. The 
trend in m pp c for polarization parallel to the c-axis is 
far less clear-cut. Band structure calculations giving the 
pressure dependence of the bare band masscj indicate a 
sublinear behaviour for its pressure dependence. How- 
ever this is a much smaller correction to linear behav- 
ior than observed here and fails to reproduce the change 
in slope at 0.5 GPa. The coincidence of 'kinks' in the 
pressure dependence of m* p obtained here and m* from 
magnetic quantum oscillation dataa is thus an indica- 
tion that the phase-phonon fitting procedure is sensitive 
to many-body effects. It is these renormalizations which 
are knowjjLto strongly influence superconductivity in this 
materialoo. 
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terms of 



dA 




din T c dlii9 



0.5 1.0 1.5 

Pressure (GPa) 



2.0 



FIG. 7: The mass mj p c for light polarized parallel to the c- 
axis, deduced from the phase-phonon fit to the optical data. 
For P < 0.5PGa, dm*°/dP ~ -2.5m e 



-58%/GPa. 
din m* pp c /dP 



dm P p/dP ~ 
For P > 0.5PGa, 
- -17%/GPa. 



dm* p c /dP 



and dlnmp P c /dP 



-0.6m c and 



Implications for the superconducting pairing 
mechanism. 



It has been shown previously^ that the usual electron- 
acoustic phonon interaction mechanism is unable to ac- 
count for the magnitude of the electron-phonon coupling 
constant or the large pressure dependence of T c in the. 
BEDT-TTF superconductors. A further refinementEZl 
includes the attractive interaction mediated by the A g 
molecular modes, with the total electron-phonon cou- 
pling constant, Atot, given by a Yamaji sum over the 
individual A g molecular modes. The energy scalg .for 
the interaction's set by the Debye frequency, OliZtEl 
Caulfield et al.a have previously inferred 9 w 40cm -1 
by fitting the effective mass dependence of T c with a lin- 
earised Eliashberg equation using an Einstein density of 
phonon states, (5(OL,High and low temperature specific 
heat measurementsE2lc2l yield values of .ranging from 
38cm -1 to 140cm -1 . Calculations of AtotlZI give values 
ranging from 0.25-0.45EM3H 

The weak-coupling BCS formulaB gives a satisfactory 
functional description of the ambient pressure T c t3 and 
describes the effective mass dependence of the supercon- 
ducting transition temperaturca. However, the fact that 
the weak-coupling BCS expression describes T c versus 
m* welfi should NOT be taken to imply that k-(BEDT- 
TTF) 2 Cu(SCN) 2 is a weak-coupling BCS superconduc- 
tor. The formula is used here as a convenient parameter- 
ization which is known to describe earlier data wellnEi 

The pressure derivative of the weak-coupling BCS for- 
mula provides a convenient parametrization of d '" ' ; " 



dP 



dP 



(A-/x*) 2 



dA d^* 
dP ~ dP 



With 



(A - A 4 *) 2 



In 



—) 

1.130 



(26) 



(27) 



m 90±50 cm -1 and using the average 

pressure-induced stiffening of the Raman active lattice 
modes of » +13 ^oGPa" 1 § for the only unknown 



is the pressure dependence of the Coulomb pseudopoten- 
tial, -Tp-. There is sufficient uncertainty in the calcu- 



lations of jti* that an accur? 
dependence cannot be made 



dP 



estimate for its pressure 
c°L However, fx* is known_to 
scale with the density of states at the Fermi energy^, 
which will decrease -ssrith the application of pressurec . 
Thus if fj>* is positives it will decrease with pressure and 
omitting ^p- from Equation ^6] will result in an overes- 
timate for the rate of suppression of T c with pressure. A 
negative value of fx* indicates that direct (non phonon- 
mediated) ialectron-electron interactions are involved in 
the pairingaca which are not measured by these experi- 
ments. Consequently ^p- is ignored in Equation p6| . 

Estimates for the pressure dependence of the electron- 
phonon coupling constant, based on the vibrations sam- 
pled and the two methods of calculation are used to cal- 
culate the right-hand side of Equation ^6|. It should be 
emphasized that the two models we are employing to es- 
timate the pressure dependence of the electron-phonon 
coupling constant lie at opposite limits of our sample's 
behaviour, the dimer charge-oscillation model being the 
limit of a localized electronic system and the phase- 
phonon model the limit of an itinerant system. A con- 
clusion common to both methods of analysis should thus 
be a robust one. 



1. T C (P) from the dimer charge- oscillation model. 

Employing the dimer charge model, Atot is assumed 
to have a pressure dependence similar to that of the 
strongly coupled C=C mode, i.e. of the order of - 
17 %GPa -1 , with an upper limit of -20%GPa -1 used in 
this calculation. This gives w -40±32 %GPa -1 , 

which is far from the experimentally observed value of 
^fp^ « -200 %GPa -1 |||. 

As can be seen from (E7 ) , logarithmically scales with 
the dependence of the pressure derivative of T c on the 
pressure derivative of the coupling constant. Thus to ob- 
tain the observed rapid fall of T c with pressure from only 
the decrease in the electron-phonon coupling constant re- 
quires to be of the order of the C=C mode frequency, 
« 1500 cm -1 . Such a value is inconsistent with the 
10 K superconducting temperature and the unconven- 
tional isotope shift observed upon carbon substitutionED. 
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Thus, the electron-phonon coupling as extracted from the 
dimer-charge-oscillation model cannot be the relevant pa- 
rameter for superconductivity in our organic supercon- 
ductor. 



2. T C (P) from the phase-phonon model. 

The pressure dependence of the electron-phonon cou- 
pling constant derived from the phase-phonon model pro- 
vides far less clear-cut results. In this case Atot is as- 
sumed to have a pressure dependence of -40±20 %GPa _1 , 
a value consistent with the phase-phonon calculation of A 
for all observed modes. This gives « -90 ^cGPa" 1 

with the asymmetric error of +150 and -77 %GPa~ 1 . The 
upper bound of this estimate is close to the experimen- 
tally observed value of sa -200 ^oGPa" 1 [§,§. 
Thus, the upper estimate of all parameters would be re- 
quired to explain the pressure dependence of T C (P) as 
arising solely due to a reduction in the electron-phonon 
coupling constant. 

3. Comparison. 

Again, it should be stressed that the dimer charge- 
oscillation model represents the limit of localized elec- 
tronic behaviour and the phase-phonon model represents 
the limit of itinerant electronic behaviour. The true prop- 
erties of k-(BEDT-TTF) 2 Cu(SCN) 2 under these experi- 
mental conditions are believed to lie somewhere between 
these extremes. Both models have similar conclusions , 
i.e. d 1 jp Tc cannot be reproduced within weak coupling 
BCS theory. The fact that the similar conclusions are 
drawn when employing either model suggests that they 
are robust. 

Use of the dimer charge-oscillation model casts consid- 
erable doubt on whether this material is a simple BCS 
superconductor because the characteristic energy of the 
pairing interaction would have to be of the order of the 
highest frequency molecular modes, a value inconsistent 
with the 10 K superconducting transition temperature. 
However this conclusion must be treated with caution 
due to the numerous assumptions involved in applying 
this model of a localized system to the coupling between 
the molecular vibrations and delocalized conduction elec- 



trons. Doubt on whether this material is a simple BCS 
superconductor is also cast by the application of the 
phase-phonon model, because the upper estimate of all 
parameters is required to explain the experimentally ob- 
served suppression of T c with pressure. 

This indicates that electron-electron interaction may 
be playing a significant role in this material's supercon- 
ducting mechanism. It is worth noting that there is 
mounting evidence in support of_this conclusion. Both 
experimentaoo and theoretically, studies predict that 
the superconducting order parameter in this material has 
an 'exotic' d-wave symmetry, a property predicted for 
a superconducthjis [-state whose pairing is mediated by 
spin fLuctuationsEj'Ell, i.e. direct electron-electron inter- 
actiontl. 



V. CONCLUSIONS. 

Comparison of high-pressure Raman scattering and in- 
frared reflectivity data provides an alternative method 
for probing the quasiparticle contributions to effective 
mass enhancement known to be intimately connected 
to superconductivity in this materiala. It also enables 
the pressure dependence of the electron-phonon coupling 
strength to be evaluated for modes observed in both 
spectra. Using the weak coupling limit of BCS the- 
ory the pressure dependence of the electron-phonon cou- 
pling constant is compared to the pressure dependence of 
the superconducting transition temperature. This casts 
doubt on whether this material is a simple BCS super- 
conductor, an indication that electron-electron interac- 
tion may be playing a significant role in this material's 
superconducting mechanism. 
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